Simulation of quantum random walks using interference of classical field 
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We suggest a theoretical scheme for the simulation of quantum random walks on a line using 
beam splitters, phase shifters and photodetectors. Our model enables us to simulate a quantum 
random walk with use of the wave nature of classical light fields. Furthermore, the proposed set-up 
allows the analysis of the effects of decoherence. The transition from a pure mean photon-number 
distribution to a classical one is studied varying the decoherence parameters. 
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I. INTRODUCTION 

Random walks are useful models for physicists to study 
statistical behaviours of nature such as Brownian mo- 
tions of free particles 0. They have also been studied 
for practical use such as algorithms in computer science 
and risk management in finance Quantum ver- 
sions of random walks have been recently studied both for 
fundamental interests and for the expectation of build- 
ing new algorithms for quantum computation Q . There 
have been several suggestions for a practical implementa- 
tion of quantum random walks, using ions in linear traps, 
optical lattices and cavity-QED Recently, pro- 

posals for the implementation of quantum random walks 
with linear optical elements have been suggested 0, E} 
and the first search algorithm using quantum random 
walks has been reported Quantum random walks 
typically show very different patterns from the Gaussian 
distributions for classical random walks, which have some 
remarkable characteristics such as an exponentially fast 
hitting time 3 . It has been pointed out that these differ- 
ences are due to the existence of quantum coherence [B| . 

In this paper, we suggest a theoretical scheme to sim- 
ulate quantum random walks on a line using the wave 
nature of classical light fields. This is related to the 
fact that the idea of quantum coherence is originally bor- 
rowed from the interference of wave mechanics shown in 
Young's double-slit experiment. In our scheme, it is also 
possible to simulate decoherence processes modeled using 
additional random phase-shifters and beam splitters with 
erratic transmittivity. This analysis is relevant under dif- 
ferent points of view. First of all it shows that, increasing 
the amount of decoherence that affects the system, the 
distribution of the random walk changes from a totally 
quantum one to a classical Gaussian distribution. This 
clarifies the role played by the interference effects in the 
dyanmics of a quantum walker and represents an ulterior 
proof of the validity of a simultation based on interfero- 
metric devices. On the other hand, studying the effects 
of possible sources of errors in our model, we can single 
out the causes of certain deviations from the ideality in 
the patterns resulting from performed experiments. 

This paper is organized as follows. In Section [HI we 
briefly review coined quantum random walks on a line 



with their characteristics and we suggest a scheme for 
the simulation of quantum random walks. We will later 
show that, with this set-up, we can simulate quantum 
random walks using wave nature of a field. This pos- 
sibility is neither always obvious nor mentioned in other 
models. It should be pointed out that a similar scheme of 
all-optical implementation has been suggested by Zhao et 
al. [?J. Their proposal is entirely based on the quantum 
coherence of a quantum superposition of two polariza- 
tion degrees of freedom. On the other hand, our scheme 
is able to take the wave nature of any input field (clas- 
sical or non-classical) to show the same interference pat- 
tern. This, with the explanation of the role of the phase 
shifters in our scheme, is shown in Section lnTl Scction llVl 
is devoted to the study of the decoherence effects in our 
proposal. We show that decoherence on the coin toss- 
ing operation and on the quantum walker motion can 
be simulated and studied by means of our system, thus 
demonstrating the role of the interference effects in this 
simulation. This investigation is useful even from a prac- 
tical point of view because it singles out and characterizes 
the effect of a class of errors that could affect a performed 
experiment. 



II. QUANTUM RANDOM WALK WITH 
LINEAR OPTICAL ELEMENTS 

In uni-dimensional coined random walks, the walker 
is restricted to move along a line with a number of dis- 
crete integer points on it. The walker is supposed to 
be a classical particle on one of the integer points. A 
coin tossing determines whether the walker moves left or 
right for each step. In the quantum version of coined ran- 
dom walks, the classical coin is replaced by a quantum 
bit whose states \L) and \R) represent the logical values 
LEFT and RIGHT. The quantum coin can be embodied 
by an internal degree of freedom of the walker itself Q . 
The walker, which is a quantum particle, moves condi- 
tioned to the result of the coin tossingoperation which 
is realized by a Hadamard transform For example, 
the transformation for one step of the particle from an 



arbitrary point X is simply 

\X, R) — -L{\X + 1,R) + \X- 1, L», 

f (1) 
^—(\X + 1,R)-\X-1,L)). 

After n steps, the state of the system is \^ n )- Differently 
from the classical walks on a line, where the position of 
the particle is monitored at every step of the process, in 
the quantum version the walker remains in a superpo- 
sition of many positions until the final measurement is 
performed. The probability for the particle being at X k 
after n steps is V n (X k ) = | (R \ (X k | *„) | 2 + 1 (L \ (X k | *„) | 2 . 
During the quantum random walk process, destructive as 
well as constructive interference may occur. The quan- 
tum correlation between two different positions on a line 
introduced at the first step may be kept by delaying the 
measurement step until the final iteration. 

The probability distribution to find the particle at a 
given position is generally dependent on the initial state 
of the system [B| and exhibits a very structured pattern. 
This allows only numerical evaluations of its variance. 
It has been shown that, roughly, the standard deviation 
oqrw grows linearly with TV and is independent from 
the initial state of the coin jy]. Thus, the walker in quan- 
tum walks explores its possible configurations faster than 
in classical walks, where the standard deviation grows 
as \N. This motivates the conjecture that algorithms 
based on quantum random walks could beat their classi- 
cal versions in terms of the time needed to solve a prob- 
lem 0. 

There have been a few suggestions for experimental im- 
plementations of quantum random walks [6j . Recently, it 
has been shown that quantum random walks can be real- 
ized using linear optical elements Q. In this scheme, po- 
larization beam splitters, half-wave plates and photode- 
tectors are used. The walker is embodied in a single- 
photon state and the entire scheme is based on the quan- 
tum coherence of two polarization states of the photon. 
This result is inspiring as a first proposal for an all-optical 
implementation of a quantum random walk, even if it re- 
quires a reliable single-photon state source, which is very 
demanding, and the apparatus is highly sensitive to vari- 
ations in the photons polarization. 

First, we propose a scheme which uses ordinary 50:50 
beam splitters, phase shifters and photodetectors. We 
formulate quantum random walks with the coin tossing 
operation embedded in the translation of the walker par- 
ticle. In our scheme, the polarization degree of freedom 
does not play a role and, thus, is not considered at all. A 
single-mode field, including a thermal field, may be used 
as an input to simulate the distribution of the quantum 
random walk. In fact, this may be apparent if we re- 
call that Young used a thermal field for his double-slit 
experiment and showed interference. 

Let us consider the experimental set-up, composed of 
50:50 beam splitters, phase shifters, and photodetectors, 
shown in Fig. ^ F° r convenience, we denote the field 
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FIG. 1: All-optical set-up for the simulation of quantum ran- 
dom walks on a line, (a) Two different kinds of operations are 
shown: Ti is an ordinary beam splitter B\(6,(f). T2 involves 
the cascade of the phase shifter Pi(n/2), of a 50:50 beam 
splitter _B2(7r/2,7r) and of the phase shifter P2(~ tt/2). (b) 
Proposed set-up, shown up to the fourth dynamic line. Apart 
the input state, all the other modes are initially prepared in 
vacuum states. 

modes propagating sidewards by s and downwards by d. 
As the beam splitters used here are polarization insensi- 
tive, these modes do not refer to polarization. Here, we 
consider a single-photon state |I) S as the initial state of 
the walker and we show that, in this case, our scheme 
gives rise to coined quantum random walks on a line. At 
the first beam splitter, £?i ($,</>), the input field is mixed 
with a field mode prepared in a vacuum state (Fig. ^a)). 
The following transformation is realized 

B 1 (9,cl))\0,l)d s = oos^\l,0) ds + e^sm^\0,l) ds , (2) 

where Bi(0,<£) = exp{6»/2(e^a|a d - e~ ili, d\a s )} is the 
beam splitter operator and a s ^ (a^ d ) are the annihila- 
tion (creation) operators for a sideward and a downward 
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field mode, respectively. We define the transformation in 
Eq. J2J) as T\. We introduce the transformation T 2 

|l,0) d , s ^-^(|l,0) + |0,l» d)S , (3) 

|0,l) < ,,.-^-^(|l,0)-|0,l» <(| . 1 (4) 

which can be realized with a 50:50 beam splitter, 
B 2 (ir/2,ir), and two phase shifters A(tt/2) = e*™^/ 2 
and P 2 (-7r/2) = e -™fad/2 as shown in Fig. [Ha). 

The scheme can simply be illustrated as recursive ap- 
plications of T 2 after the initial transformation T\, as 
shown in Fig. db). A dynamic line 7] is represented 
by a row of aligned optical elements (or photodetectors), 
labelled j in Fig. CJb). On the other hand, a node is 
given by a point represented by A: on a dynamic line. For 
example, the detector £>_ 2 is on the fourth dynamic line 
and occupies the node k = — 2. If a photon is incident 
downward (sideward) on a dynamic line j and node k, 
we represent its state as \k, d)j (\k,s)j). The transition 
from a dynamic line j to j + 1 by means of the operation 
T 2 is synthesized by 

\k,d)j^ ^(\k+l,d) + \k-l,s)) j+1 , 

\ ( 5 ) 
\k,s) J ^—(\k + l.d)-\k-l, S )) 3+1 . 

We notice that Eqs. JSJ are equivalent to Eqs. (|TJ- Thus, 
the actions of T± and T 2 on a single-photon state exactly 
corresponds to a coined quantum random walk. Any 
initial coin state, up to an irrelevant global phase, can 
be prepared changing 9 and ip in T\. If 6 = vr/2 and 
4> = —tt/2, we get the symmetric probability distribution 
that corresponds to the initial coin state (\R)+i\L))/y/2 
in a coined quantum walk [6|. In our model, the differ- 
ence between quantum and classical walks from a certain 
step is due to the interference of the walker's paths on 
the T 2 processes 



does not lead to any entanglement between the output 
modes 13]. Assuming 9 — n/2 and <j) = — 7T /2 for the Xi 
process, we can calculate the distribution of the average 
photon-number as a function of the position k on the 
chosen final dynamic line. For example, for N = 4 steps, 
we find the final state 




(6) 

with \a, s) ■ (\a,d)j) that indicates a coherent state inci- 
dent sideward (downward) on a dynamic line j and node 
k. The average photon-number Af p (N,k) for node k is 
Af p (4, k) = M(A, k)Af in (\a)), with 

M(4,±4) = -L M(4,±2) = ~ M(4,0) = ~. (7) 

Here, A/i„(|a)) = |a| 2 is the average photon number 
for the input state \a) and M(N,k) is the normalized 
photon-number distribution at step N and node k. It 
characterizes the output photon-number distribution at 
the detectors. We find that the distribution A4(A,k) for 
an input coherent state is the same as the one for the sin- 
gle photon input : 6], i.e., the two different inputs result in 
the same photon-number distribution. The average pho- 
ton numbers for steps, N = 4, 5, 6 are shown in Fig. [21 
The deviations of a quantum walk from its classical coun- 
terpart appears from the fourth step. This is due to the 
particular values of the parameters in the transformation 
T\\ 9 — it/2 and <fi — —ir/2. Since a coherent state input 
results in the same quantum random walk pattern of the 
single photon case for all the steps we have considered, 
we conjecture that the quantum walk pattern results even 
any initial state for a general number of steps N. In what 
follows, we prove the validity of this conjecture. 

B. General case 



III. ANALYSIS WITH DIFFERENT STATES OF 
THE WALKER 

In this Section we show that the scheme suggested in 
Fig.ncxhibits the same interference pattern at the detec- 
tors regardless of the nature of the input state. We first 
address the case of an input coherent state and, then, we 
extend the analysis to any field. 

A. Coherent states 

A coherent state \a) (a G C) is generally assumed to 
be the best description of the state of a laser beam. We 
consider \a) as the input state of the walker. The action 
of the beam splitter operator on two input coherent states 



With the proposed set-up, the quantum walk process 
can be represented as 

\$ N ) = UTU=N):.UTU=l)f 1U=0) \^o) = U$ w \®o), 

where |$q) is the input state, N is the number of steps, 
and Ut is an appropriate unitary transformation for each 
step. For a coherent state, the previous result can be 
summarized as 

jjN 

l$o) = \a) \xia)i\x2a) 2 --\XNa)2N = |$jv)< (8) 

Eq. is an explicit example. The average photon num- 
ber for mode r (0 < r < 2N) is n r = IxrPM 2 = 
|Xr| 2 Mn(|a)), w hh r — corresponding to the mode 
incident on the detector that occupies j = N, k = —N. 
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FIG. 2: Average photon- number distribution for an input co- 
herent state | a = 1), as a function of the position along the 
final dynamic line. Three different cases are considered: the 
solid-line curve is relative to a number of steps N = 4; the 
dashed-line represents N = 5 while the dot-dashed one is for 
TV = 6. The plots match perfectly the graphs expected for a 
coined quantum walk on a line. In the general case of a / 1, 
N p has to be normalized with respect to \a\ 2 . 

It is easy to show that the average photon number for 
the fc-th node and j-th step is given by Af p (j,k) = 
Uj-k +rij_fc + i = (\xj-k\ 2 + \X]-k+i\ 2 Wtot(\a)), where 
Xo = X2N+1 = 0. This result also means that 

M(j,k) = \ Xj - k \ 2 + \Xj-k + i\ 2 . (9) 

Note that Xr does not depend on the amplitude of the 
initial state but only on the structure of Uq W . 

The initial state density operator in P representation 
can be generally written as 0, [lj] 

Po = J d 2 aP(a)\a)(a\ 1 (10) 

where P(a) is the P representation of the initial state po- 
Provided that P(a) is a sufficiently singular generalized 
function, such a representation exists for any given oper- 
ator po 13- After N steps, the density operator evolves 
as: 

Pn = Uq W p Uq W 

= J d 2 aP{a)\x\a) i{xia\ ® .. ® \x2Na) 2 N(X2Na\ 

(11) 

where Eqs. (JSJ) and (|10(l have been used. The P repre- 
sentation is particularly appropriate for our aim to find 
the average photon-number distribution since it can be 
shown that the moments of the P representation give 
the expectation values of normally-ordered products of 
bosonic operators [jllT^. 

The marginal density matrix for mode r is simply ob- 
tained as 

p r = / d 2 aP(a)\x r a) r (Xra\. (12) 



The average photon number for the r-th mode is 

n r = Tr r [p r a^a] = \x r \ 2 J d 2 aP(a) \a\ 2 = |Xr| 2 W m (po), 

and the average photon number for the j-th step and 
fe-th node is Af p (j,k) = M(j, fc)Mot(po) = (\Xj-k\ 2 + 
\Xj~k+i\ 2 )Ntot(po), from which Eq. © is found to hold 
for the case of any input field. The interference pattern 
determined by M.(j, k) does not depend on the initial in- 
put state. For a given set of beam splitters and phase 
shifters, any input state will result in the same interfer- 
ence pattern. Only an overall factor will be changed, ac- 
cording to the total average photon-number of the initial 
state. For a classical light, in a pictorial way, the result 
is nothing but quantum random walks with many walk- 
ers simulated by interference between fields. For a weak 
field, the quantum random walks with a single walker 
can be probabilistically performed. For example, given 
a coherent state with a = 1, a single photon is detected 
with 37% of the probability. 

A problem of the approach employing dynamic lines 
for quantum random walks is that the required number 
of resources (in terms of the number of optical elements 
required for a chosen number of steps and of the field 
modes involved) grows quadratically with the number 
of steps. This imposes serious limitations to the scal- 
ability of such a proposal and affects the efficiency of a 
simulation based on an interferometric device. In the 
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FIG. 3: Alternative set-up for quantum random walk on a 
line. In this scheme, the number of required resources scales 
linearly with the nymber of steps N. Two rows of Acousto- 
Optic Modulators (AOM) direct the incoming beams of light 
to the perfect mirrors M or to the detectors row. This set-up 
is conceptually equivalent to that sketched in Fig. 0^b). 

alternative proposal in Fig. |3 this problem is bypassed 
measuring all the even positions by the upper row of de- 
tectors, while the odd ones are detected by the lower 
row. Acousto-Optic Modulators (AOMs) are used to 
guide a beam toward a mirror for further steps or toward 
a detector for the measurement. When the AOMs in the 
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top row have to deflect the light beams toward the detec- 
tors, those in the bottom row should not be active. The 
beam splitters and phase shifters in Fig.Qand in Fig. 
are the same. The number of required resources, in this 
latter scheme, increases only linearly with the number of 
steps [lflj . 

There are many difficulties, for a practical implemen- 
tation of the schemes we propose, that have to be taken 
in consideration. For example, being a multimode in- 
terferometric apparatus, the proposed set-up could be 
affected by the misalignment of the involved optical ele- 
ments. Furthermore, we need 2N modes for N steps of 
the walk process, that makes the controllability of the 
system very difficult. Nevertheless, even if these prob- 
lems render the proposed set-up challenging under an 
experimental point of view, our proposal has to be seen 
as a thought experiment useful for the investigation of 
the physics that is behind the appearance of the charac- 
teristic probability pattern of a quantum random walk. 



IV. SIMULATION OF DECOHERENCE IN 
QUANTUM RANDOM WALKS 

To better understand how interference effects are at 
the basis of a quantum walk process we study the effect 
that a certain class of errors have on the performance of 
the set-up we propose. A decoherence mechanism is po- 
tentially able to wash out the interference pattern, thus 
erasing the speed-up characteristic of a quantum walk 
and restoring some aspects of the classical diffusion pro- 
cess. In this Section we study two different models for 
decoherence in our set-up. We show the transition of 
the dynamics of the walker from the pure quantum to 
the classical case. This analysis is in part motivated by 
the attention that has been recently payed to the way in 
which the quantum walk pattern is modified by imperfect 
coin tossings or walker translations, both for quantum 
walk on a line and higher dimensions 0, ll(| . A remark- 
able result, shown by Kendon and Tregenna in |10|. is 
that small amounts of decoherence, rather than render 
the process useless for the purposes of quantum informa- 
tion, amazingly increase the capability of the system to 
explore its possible configurations. This gives a proba- 
bility distribution to find the walker in a certain position 
that spreads faster than in pure dynamics. Our study 
is able to highlight even this aspect of the dynamics of 
the walker. On the other hand, studying the effects of 
possible sources of decoherence is worth under a practi- 
cal point of view. The characterization of some relevant 
sources of errors, in the proposed set-up, will make us 
understand why the pattern resulting from a performed 
experiment devoted to the realization of a quantum walk 
process could deviate from the ideality. 

We have considered ulterior phase shift operations per- 
formed just before and after each Ti transformation. The 
shift in these additional operations is randomly chosen 
from a Gaussian distribution. In what follows, we show 



how the mean photon-number distribution changes its 
shape (from a classical Gaussian pattern to an approxi- 
mately flat distribution then to a quantum distribution) 
as the amount of randomness in the additional phase 
shifts is reduced. 

If I is a number randomly taken from a Gaussian distri- 
bution centred at 1 with an adjustable standard deviation 
a pp , we shift the phase of each field mode in Fig. ^b) 
by an amount equal to 2tt If the phase shift is equal 
to 2ir, the additional phase shifters are ineffective and a 
quantum walk pattern is recovered. On the other hand, 
if the amount of shifts deviates from this neutral value, 
they affect the interferences responsible for the quantum 
walk and some deviations have to be expected and the 
average over a large number of trials results in a classical 
distribution. 
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FIG. 4: Average photon-number distribution vs position for 
an input coherent state \a — 1} and for 200 steps. Different 
cases are considered: the bell-like curve represents the case of 
an introduced randomness I taken from a Gaussian distribu- 
tion centred at 1 and with a vv = 0.25. The curve evidently 
resembles the expected Gaussian distribution. The solid line 
shows the results for a pp = 0.0125. It is compared with a uni- 
form distribution between — 200/\/2 and 200/\/2, fictitiously 
extended to improve visibility. Finally, the dashed curve rep- 
resents the pure quantum case corresponding to / chosen from 
a Dirac delta function 5(1 — 1). Each point in the simulated 
curves is averaged over 50 different trials. 



In Fig. QJ the shown distributions are the results of an 
average over 50 different trials: in each one of them, and 
for each step in a single trial, a different random value 
for I is considered and the mean photon-number at the 
various locations on the final dynamic line is calculated, 
averaging over the outcome for each trial. 

If, now, (T pp is reduced (in Fig. a pp — 0.0125), the 
phase shifts vary over a small range of values around 
2ir. The dynamic evolution of the system is affected 
in such a way that no classical signature is evident in 
the mean photon-number distribution. A highly non- 
classical pattern is found and some deviations from the 
pure quantum random walk case are evident. The distri- 
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bution is relatively flat over a region that is wider than 
the pure quantum case. This result is in good agreement 
with the analysis performed in [Tol ] for a small amount 
of decoherence. In our case, the limited randomness im- 
posed to the evolution of the photonic walker simulates 
the effect of a decoherent coin tossing. The remarkable 
feature in this analysis is that we have used just classi- 
cal resources (linear optics elements and input coherent 
states). Nonetheless, we still simulate the relevant fea- 
tures of the transition from a pure quantum evolution to 
the classical spread due to a large superimposed random- 
ness. 



In [X5J , the effect of phase randomness in a general in- 
tcrfcromctric device has been investigated. In particular, 
if the device can be thought as the iterative applications 
of some basic units, each one affected by a fixed random- 
ness ^5|, then Anderson localization can be obtained. 
Indeed, when fixed randomness is considered, a connec- 
tion to the th eory of the band-diagonal transfer matrix 
(examined in [Tg) can be established. It is this kind 
of dynamic evolution that leads to localization of the 
walker. Physically, the model described in this case is 
near to the repeated passages of a beam of light through 
a dielectric layer placed inside an electro-magnetic cavity, 
as described in [17|. In our model, however, no localiza- 
tion effect is achieved since different values for the phase 
shifts at each step are taken. In this respect, our case is 
far from a band-diagonal evolution. These qualitative ar- 
guments are resumed in Fig. [5[ where the transition from 
a flat distribution (obtained for a small decoherence pa- 
rameter (jpp) to the classical one (relative to a strongly 
randomized quantum walk) is reported. To compare our 
results to those in 0] and to show that no dynamic lo- 
calization is here achieved, we present plots for the aver- 
age photon-number distributions in lin-log and in lin-lin 
scale. 



Following the same lines depicted above, we can in- 
vestigate about errors due to the uncertainty in the 
beam splitters transmittivities. We consider imperfect 
beam splitters whose transmittivities randomly fluctuate 
around 50% according to a Gaussian distribution with 
standard deviation er^. Computing the normalized av- 
erage photon-number distribution for an input coherent 
state, we find a narrow range of values for er^ within 
which a flat distribution is achieved. Outside this range, 
the distribution rapidly converges toward a classical one. 

To give a picture of the combined effect of the two deco- 
herence processes, we include random phase shifters be- 
tween two subsequent T2 operations and random fluctua- 
tions in the transmittivity of the beam splitters. In Fig.[f)| 
we show the distribution that corresponds to a pv = 0.005 
and <7&s = 0.07. We can see that the mean-photon distri- 
bution has been very much flattened. Of course, as a pp 
and (Tbs grow, the curve will become Gaussian. 
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FIG. 5: Transition from weak to strong randomization in the 
model for decoherence in the coin tossing for an input coherent 
state \a — 1). From top to bottom, a pp , is increased. We have 
considered a pp — 0.013 (a), a pp — 0.13 (b) and a pp — 0.25 (c). 
The figures in the right show the same distributions presented 
in the left but in lin-log scale, with which the investigation 
of the appearance of localization effects is easier. The mean 
photon-number distribution smoothly changes from a sharp- 
squared distribution to a concave curve that is typical of a 
classical distribution flEf. 



V. REMARKS AND DISCUSSION 

As we discussed, the realization of the model we pro- 
pose is not trivial as we pays the price represented by 
the use of 2N field modes to replace the quantum walker 
(that belongs to a Hilbert space of dimension N). Thus, 
the addition of a component, in our set-up, increases the 
difficulty of alignment. However, what we want to stress 
in this paper is the possibility of simulating quantum ran- 
dom walks using the wave nature of a classical field. We 
have shown that this study has been possible using our 
thought experimental set-up. Even though this results 
could be surprising, the possibility of such a simulation 
may be a natural result if we consider that quantum co- 
herence and quantum interference are concept originally 
borrowed from wave mechanics. This possibility has been 
formally proved using standard tools of quantum optics. 

Furthermore, we have simulated some decoherence 
mechanisms on the quantum random walk by means of 
linear optical devices and input coherent states. We have 
observed how the average photon-number distributions 
are modified when controlled randomness is introduced 
in the system via additional phase shifters and imperfect 
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the coin to be in the \R) or in the \L) state As we 
P have seen in Section^ properly choosing the parameters 

of the optical elements in Ti, T2, our proposal is able to 




FIG. 6: Average photon-number distribution for an input 
state \a = 1) considering both the models of decoherence. 
The number of steps considered is 200 and each point is av- 
eraged over 50 different trials. We have taken a pp = 0.005. 
On the other hand, we have taken = f |m| in T2, with a 
random number, m, extracted from a Gaussian distribution 
centered at 1 and having standard deviation crbs = 0.07. 

beam splitters. This analysis is useful both theoretically 
(clarifiying the role of the coherent effects in the simu- 
lation) and practically because it characterizes the influ- 
ences of possible sources of errors affecting the results of 
a performed experiment. 

Finally, we want to mention here that it is in principle 
possible to extend our scheme to quantum random walks 
on a circle of N points, as shown schematically in Fig. Q 
One can adapt the concept of dynamic line to that of 
dynamic circles: the walker transits from circle to circle 
(each having a non decreasing number of sites on it) in 
a way completely similar to that described in Section ITU 
The number of required dynamic circles is equal to N . 
Each site on a given circle is occupied by a basic oper- 
ation: T\ occupies the unique site on the first dynamic 
circle, all the other sites in the following circles (labeled 
as j = 1, .,,N in Fig. 01 being occupied by T 2 . After each 
Ti operation, the beams are directed, by means of some 
mirrors, toward the proper site on the next dynamic cir- 
cle, as shown in Fig. for the transition from the j = 
to the j = 1 circle. At the final dynamic circle, the mean- 
photon number distribution at the sites is revealed by an 
array of detectors. Basically, this implementation is still 
based on the simulation of a quantum walk on a line and 
it is, thus, obvious that it will simulate quantum walks 
on a circle with classical fields. 

Under certain circumstances, our approach can be use- 
ful in order to simulate quantum walks on a hypercube 
of dimension 3. This higher-dimensional quantum walk 
can, indeed, be reduced to a biased quantum walk on a 
line with properly chosen, asymmetrical, probabilities for 




FIG. 7: An implementation of a quantum random walk on a 
circle using dynamic circles. 



realize quantum random walk on a line with any biased 
coin. We, thus, expect the possibility to simulate quan- 
tum walks on a three-dimensional hypercube by means of 
interference of classical light. However, the extension of 
these results to general graphs and hypercubes of higher 
dimensions (as well as an analysis of the efficiency of such 
a simulation) is much more difficult and goes beyond the 
purposes of this work. It is, however, worth stressing 
that the possibility of classical simulations of quantum 
random walk on a line does not necessarily imply their 
usefulness for a practical quantum algorithm. There re- 
main important open questions about the gain, in terms 
of speed-up of quantum computation, that can be ob- 
tained from the classical simulation of quantum walks. 

Note added - Knight et al. also pointed out the 
possibility of simulation of quantum random walks using 
classical fields but using a totally different set-up [l8| . 
This appeared one day before the present manuscript of 
ours was uploaded into the Los Alamos archive. 
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